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I.  INTRODUCTION 

Materials  such  as  propellants  and  explosives,  must  by  their 
nature  contain  tightly  bonded  groups  of  atoms  which  retain  their 
molecular  character  until  a  sufficient  stimulus  Is  applied  to  cause 
dissociation.  For  the  material  to  be  "energetic",  the  dissociation  is 
accompanied  by  an  energy  release  which  in  turn  triggers  further  dis¬ 
sociation  leading  to  Initiation  or  ignition.  The  desired  macroscopic 
behavior  is  ultimately  controlled  by  the  Interatomic  forces.  Only 
when  it  is  possible  to  quantitatively  describe  the  forces  will  it  be 
possible,  to  predict  whether  a  given  molecular  structure  will  support 
an  explosive  reaction.  The  present  inability  to  make  such  predictions 
represents  one  of  the  fundamental  limitations  to  the  advance  of  energe¬ 
tic  materials  science  and  technology. 

At  the  1970  and  1972  Army  Science  Conferences  the  authors  of 
the  present  paper  demonstrated  that  some  of  the  dynamic,  molecular  pro¬ 
perties  of  materials  of  the  complexity  of  explosives  and  propellants 
could  be  elucidated  if  a  comprehensive  experimental  and  theoretical 
approach  were  employed  (1,2).  The  present  paper  represents  Implementa¬ 
tion  and  further  refinement  of  the  capabilities  required  to  character¬ 
ize  the  fundamental  processes  of  Interest  in  energetic  materials.  In 
particular,  attention  is  focused  on  the  correlation  of  interatomic 
forces  and  physical  and  chemical  instability  as  manifested  by  changes 
in  the  dynamics  of  molecular  groups  in  ammonium  nitrate,  ammonium  per¬ 
chlorate  and  the  metal-azide  salts. 
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Ammonium  nitrate  (NH^NO-j)  Is  of  considerable  military  Interest 
because  of  its  explosive  properties.  furthermore  it  Is  cheap  and 
available  in  enormous  quantities  so  that  its  more  wide-spread  military 
use  could,  potentially,  result  in  cost  savings  of  many  millions  of  dol¬ 
lars.  The  limited  success  of  conventional  approaches  in  providing  a 
solution  to  problems  associated  with  its  use  has  lead  to  their 
examination  with  the  techniques  of  the  present  program. 

Ammonium  perchlorate  (NH4CIO4)  is  a  compound  used  extensively 
in  propellant  formulations  and  mav  well  be  the  most  thoroughly  studied 
of  all  energetic  materials.  This  because  of  its  military  usefulness, 
hut  also  because  aspects  of  its  behavior  have  general  scientific  impor¬ 
tance.  Despite  the  extensive  study  no  coherent  picture  of  the  inter¬ 
relationship  between  its  molecular  properties  and  macroscopic  behavior 
had  yet  emerged.  The  results  of  the  present  investigation  go  far  to 
rectify  this  failure. 

The  metal-azide  salts  form  a  family  of  compounds  exhibiting 
chemical  properties  ranging  from  stable  to  highly  unstable.  For  this 
reason  comparison  of  the  properties  of  different  metal-azide  salts  and 
related  compounds  should  provide  some  insight  into  the  critical  para¬ 
meters  involved  in  dissociation  and  its  propagation. 

Results  for  the  compounds  studied  are  presented  in  Section  Ill. 
These  are  preceded  in  Section  II  by  a  discussion  of  experimental  tech¬ 
niques  including  that  of  quasi-elastic  neutron  scattering  which  has  not 
previously  been  employed  for  energetic  materials.  A  brief  introduction 
to  the  theory  of  lattice  or  molecular  dynamics  is  also  presented  in 
Section  II. 


II.  THEORY  AND  EXPERIMENTAL 


Lattice  dynamics  is  the  term  generally  applied  to  the  quantitative 
description  of  atomic  and  molecular  motions  in  crystals  and  the  forces 
of  interaction  which  produce  the  motions.  The  usual  method  of  expres¬ 
sing  the  equations  of  motions  of  the  particles  in  a  crystal  is  in  terms 
of  cartesian  displacement  coordinates  X11.  With  this  coordinate  system, 
the  potential  energy  of  the  crystal  in  the  adiabatic  and  harmonic 
approx'mation  is  written  as  (3) : 
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where  a  numbers  the  S  particles  in  the  nth  primitive  cell.  The  pro¬ 
perties  of  the  force  constant  matrix  I  are  dependent  upon  the  type  ot 
force  field  used  to  describe  the  interaction  of  the  various  particles 
and  must  be  consistent  with  the  svmmetrv  of  the  crystal.  The  equations 
of  motion  to  be  solved  are 
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where  uij(q)  and  Kq  are,  respectively,  the  vibrational  frequency  and 
quasi  momentum  of  the  jth  mode  phonon,  n  ■  n  -  n',  and  Mq,  is  the  mass  of 
the  ath  atom.  A  "phonon"  is  the  quantized  plane  wave  of  wavelength  \  m 
2tt / q  and  energy  E  *  Koj j ,  in  terms  of  which  the  vibration  is  described. 

The  quantity  LajC<J)»  called  the  polarization  vector,  describes  the 
magnitude  and  pnase  of  the  motion  of  the  ath  particle  when  the  phonon 
(j,q)  is  excited,  ’’’he  determination  of  to j  (q)  and  £,tj(q)  characterizes 
a  vibration,  and  the  functional  dependence  of  these  quantities  is 
determined  by  crystal  symmetry  and  the  force  field  between  particles. 

The  measurement  of  these  two  functions  is  the  aim  of  experimental 
lattice  dynamics. 

The  point  to  be  emphasized  with  regard  to  the  conventional 
formulation  of  the  description  of  lattice  dynamics,  outlined  above,  is 
that  the  atoms  in  reao.  crystals  d£  not  experience  a  harmonic  potential. 
The  degree  of  anharmonicity  is  different  for  each  type  of  atom  and  for 
each  direction  in  which  the  atom  moves.  The  importance  of  anharmonicity 
in  the  potential  can  be  better  understood  when  it  is  realized,  for 
example,  that  in  a  purely  harmonic  crystal  there  is  no  thermal  expansion, 
phonons  do  not  interact,  and  phase  transitions  cannot  occur.  Of  course, 
all  of  these  occur  in  real  crystals. 

The  degree  of  anharmonicity  in  the  interatomic  potential  is 
best  determined  by  the  study  of  the  changes  in  vibrational  mode  frequen¬ 
cies  as  a  function  of  temperature  and  pressure.  Anharmonicity  affects 
the  frequencies  in  two  different,  ways.  One  of  these  is  a  purely  volume 
dependent  shift.  The  other  shift  -  called  the  "self-energy"  -  arises 
from  the  direct  interaction  of  the  vibrational  modes,  and  is  both  volume 
and  temperature  dependent.  Experimentally,  these  anharmonic  contribu¬ 
tions  can  be  separated  only  through  a  combination  of  spectroscopic 
studies  performed  at  various  temperatures  and  pressures,  coupled  with 
pressure-volume  and  temperature-volume  measurements. 

The  current  theory  (4)  which  describes  anharmonic  effects, 
although  formally  straightforward,  has  been  limited  in  application 
to  a  very  few  simple  solids;  details  will  not  be  presented  here.  In 
cases  of  large  anharmonicities,  the  theory  has  not  been  fully  developed. 

°ne  aspect  of  the  general  case  of  large  anharmonic  effects 
which  has  recently  gained  increasing  importance,  particularly  in  the 
case  of  ferroelectric  materials,  Is  that  of  displacive  phase  transitions. 
It  has  been  found  that  a  ferroelectric  to  paraelectric  phase  transition 
is  often  precipitated  by  a  particular  translational  mode  going  "soft" 
(i.e.  wj  (4) ^>o  as  T-*TC)  as  a  consequence  of  increasingly  large  amplitude 
vibration  resulting  from  anharmonicity. 
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Rotatory  motions  of  molecular  groups  (l.e.  motions  of  atoms 
In  a  molecule  about  the  center  of  mass  of  the  molecule  itself)  are  also 
expected  to  manifest  the  effects  of  anhannonlc  interactions.  A  direct 
consequence  of  the  anharmonic  potential  in  real  crystals  is  the  fact 
that  for  a  molecular  group,  a  rotatory  motion  of  sufficiently  large 
amplitude  can  result  in  a  "permanent"  reorientation  of  atoms  in  the 
molecule  to  equilibrium  sites  vacated  by  atoms  of  the  same  molecule. 

This  is  illustrated  in  Fig.  1  for  an  NH^  ion.  Only  recently  has  this 
type  of  anharmonic  effect  come  under  quantitative  examination. 

Interatomic  forces  cannot  be  measured  directly  but  are  Inferred 
from  measurement  of  the  to j  (<f)  and  tajCf).  (1,2)  The  principal  experi¬ 
mental  techniques  used  for  lattice-dynamical  studies  are  Raman  and  infra¬ 
red  spectroscopy,  and  coherent  neutron  inelastic  scattering  (CNIS). 

Raman  and  infra-red  spectroscopy  provide  information  only  for  u>j  (q“o) , 
whereas  CNIS  is  not  restricted  in  q.  However,  not  all  vibrational  modes 
are  observed  by  each  experimental  technique,  but  are  subject  to  "selec¬ 
tion  rules".  (2)  In  addition,  neutron  or  x-ray  diffraction,  by  which 
the  crystal  structure  is  obtained,  are  necessary  pre-requisites  to 
lattice-dynamical  studies.  In  solids  of  the  complexity  of  energetic 
materials,  all  experimental  probes  provide  Important  information  for 
determination  of  the  interatomic  forces. 

A  relatively  new  experimental  technique  -  that  of  quasi-elastic 
neutron  scattering  (QNS)  -  has  been  employed  for  energetic  materials  for 
the  first  time.  This  type  of  probe  gives  direct  information  for  the 
rotatory  reorientations  of  molecular  groups  not  readily  obtainable  by 
other  means.  Consequently,  QNS  provides  a  means  of  characterizing  the 
appropriate  interatomic  potentials  for  this  special  type  of  dynamical 
behavior,- which  can  range  from  free  rotation  to  small  amplitude  rotatory 
vibration  with  no  reorientation. 

The  QNS  technique  is  sensitive  to  this  entire  range  of  behavior; 
however,  it  is  the  intermediate  region  consisting  of  rotatory  vibrations 
for  a  period  of  time,  t,  long  compared  to  a  vibrational  period,  followed 
by  a  reorientation  which  is  of  greatest  Interest,  t  is  generally 
referred  to  as  the  residence  time.  ^ 

In  general,  the  QNS  cross-section,  o(Q,u>)  for  an  atom,  is  given 

by  (3) 


o(Q.u)'v  JdrJ  dt  exp[i(Q*r-wt)]^Pi(t)6(r-r1) 

— ► 

where  Q  and  Ku  are,  respectively,  the  momentum  and  energy  transferred 
by  the  neutron  in  the  scattering  process,  r^  denotes  the  position  of  the 
scattering  atom,  and  Pj(t)  is  the  probability  of  finding  the  atom  at 
site  i  at  time  t  given  that  it  was  at  rc-o  at  t-o.  In  the  temporal 
region  of  interest,  one  can  assume  that  the  rotational  jumps  are  instan¬ 
taneous  so  that  the  particle  la  always  at  one  of  the  equilibrium  sites. 
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The  Pi(t)  are  then  obtained  in  a  relatively  straightforward  manner. 
Skold  has  done  this  for  a  tetrahedral  molecule  with  all  four  sites 
equivalent,  l.e.  equal  x's.  (5) 

We  have  extended  the  theory  to  the  case  more  appropriate  to 
NH4CIO4:  t04ti“T2“T3.  This  residence  time  distribution  yields  a  QNS 
cross-section  (6) 


a(Q,a,)^ir[l+3j0]  6(ui)  +  [l-JollU^)  +  L(a2)  +  L(a3)] 


where  jQ  -  slnQR/QR,  <S (u>)  is  the  Dirac  delta  function, 


L(ctj)  -  aj/(a2  +  w2) 
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and  R  Is  the  jump  distance;  in  ammonium,  R  ■  1.78  A. 

In  practice,  the  measured  o(Q*oi)  consists  of  a  resolution  broadened 
6(u)  contribution  superimposed  on  the  Lorentzlan  components.  The  full- 
width-at-half-maximum  (FWHM)  of  the  composite  distribution  vs.  Q  gives 
a  crude  estimate  of  the  character  of  the  reorientation  and  the  ts;  a 
detailed  analysis  of  the  line-shapes  at  each  Q  Is  a  more  sensitive  test 
of  the  model.  Both  approaches  are  used  In  the  next  Section. 


III.  RESULTS 


A.  NH4NO3 

Ammonium  nitrate  crystallizes  In  as  many  as  five  different 
polymorphic  forms  at  atmospheric  pressure,  and  It  is  their  existence 
that  leads  to  serious  technological  problems  for  the  utilization  of 
NH4NO3  in  military  hardware.  For  example,  the  phase  transition  between 
phase  IV  (-18#C  +-*■  32.3*C)  and  phase  III  (32.3°C  *-►  8A.2*C)  occurs  at 
90®F  which  is  a  temperature  through  which  munitions  In  storage,  during 
loading,  or  in  the  field  are  cycled  by  climatic  and  manufacturing 
conditions.  One  practical  problem  which  arises  is  that  the  temperature 
coefficient  for  expansion  of  phase  III  Is  much  larger  than  that  of 
phase  IV  so  extrusion  of  NH4NO3  through  casings  occurs,  which.  In  turn, 
alters  performance.  The  Interest  here  Is  to  determine  the  fundamental 
properties  of  NH4NO3  in  Its  various  phases,  which  will  suggest  methods 
by  which,  for  example,  the  above  problem  can  be  eliminated. 
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At  the  1972  Army  Science  Conference,  the  first  complete  struc¬ 
tural  study  of  NH4NO3  (IV)  was  reported.  (2)  Comparison  of  the  detailed 
structure  of  phase  IV  and  the  partial  determinations  of  the  p  \ase  II 
and  IV  structures  revealed  that  the  occurrence  of  phase  III  was,  in 
fact,  anomalous.  That  is,  a  structural  transition  IV  ■*-*■  II  involved 
much  less  rearrangement  and  less  energetic  processes  than  IV  *-*■  III  or 
II  ■*-*  III.  From  the  microscopic  point  of  view,  a  catalyst  -  known  to 
be  surface  or  occluded  water  (2,7)  -  was  required  for  phase  III  to  occur 
To  attempt  to  shed  further  light  on  this  problem  QNS  was  employed,  for 
the  first  time,  to  examine  the  rotational  dynamics  of  NH4  in  ammonium 
nitrate.  (6) 

Phases  IV,  III,  II  (84.2  +-*■  125#C)  and  I  (125°  <-»  MP)  were 
studied  after  mild  drying  at  100°C  ("wet")  and  after  drying  at  135°C 
("dry").  Approximately  IX  additional  H2O,  by  weight,  was  removed  after 
redrying  at  the  higher  temperature.  Powder  neutron  diffraction  measure¬ 
ments  were  performed  as  a  function  of  temperature  to  confirm  earlier 
phase  transition  findings.  QNS  studies  were  made  at  temperatures  of  24° 
45°,  75°,  95®,  120°  and  130®C,  wet  and  dry.  Line  shape  analyses  of  the 
individual  spectra  have  not  been  completed;  however,  FWHM  vs  Q  plots 
are  shown  in  Fig.  2  for  phases  IV  and  IT,  wet  and  dry.  (6)  The  follow¬ 
ing  conclusions  can  be  drawn: 

1.  The  powder  diffraction  measurements  confirm  the  occurrence 
of  phases  IV,  III,  II  and  I  in  wet  samples.  Furthermore,  drying  at 
135®C  removes  sufficient  water  to  eliminate  phase  III  so  that  the 

IV  ■*-*  II  transition  occurs  directly,  at  ^50°C. 

2.  The  presence  of  a  small  amount  of  H2O  clearly  affects  the 
NH4  dynamics  in  phase  II.  The  nature  of  the  change  in  the  FWHM  vs  Q 
curves  for  phase  II  indicates  shorter  residence  times  and  weaker  (NH4)- 
(NO3)  bonds  in  the  wet  sample  than  in  the  dry.  It  is  believed  that 
phase  IV  wet  and  dry  wo>:ld  exhibit  similar  differences  except  that  they 
are  masked  in  our  measurements  by  statistical  uncertainties  and  insuffi¬ 
ciently  fine  resolution. 

3.  Surface  (or  occluded)  water  in  NH4NO3  causes  a  cooperative 
weakening  of  N-H*’*0  hydrogen  bonds  throughout  the  bulk  of  the  material. 
This  weakening  of  H-bonds  causes  phase  IV  to  become  unstable  at  a 
temperature  where  phase  II  is  not  yet  energetically  favorable.  This 
results  in  the  formation  of  the  anomalous  phase  III. 

B.  NH4CIO4 

Ammonium  perchlorate  is,  as  indicated  above,  a  most  thoroughly 
studied  material;  nevertheless,  the  macroscopic  properties  are  still 
far  from  understood.  In  particular,  the  peculiar  decomposition  behavior 
(8)  in  the  neighborhood  of  the  orthorhombic  to  cubic  phase  transition 
(513°K)  is  puzzling  but  clearly  indicates  that  the  crystal  structure 
and  the  stability  of  intramolecular  bonds  are  closely  inter-related. 

It  is  toward  elucidating  the  microscopic  properties  which  cause  this 
oehavior  that  our  studies  of  NH4CIO4  are  ultimately  directed.  At  the 
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same  time,  the  low  temperature  properties  of  NH4CIO4,  although  the 
subject  of  extensive  scientific  study  for  several  years,  have  also 
remained  uncertain.  (9) 

To  clarify  this  long-standing  problem,  and  to  establish  a 
framework  within  which  to  examine  the  high-temperature  behavior,  the 
microscopic  properties  have  been  examined  in  considerable  detail  in  the 
10#-300°K  temperature  range  by  means  of  neutron  diffraction,  QNS  and 
Raman  scattering. 

Previous  experimental  work  (i.e.  NMR,  Cv,  IR,  Raman  and  neutron 
transmission)  has  been  analyzed  to  characterize  NH^  motions.  (9)  Based 
on  assumptions  that  the  reorientations  occur  entirely  with  one  hydrogen 
fixed  or  with  all  four  hydrogens  bound  equally  (i.e.  equal  t's),  these 
studies  have  yielded  potential  barriers  ranging  from  0.2  to  2.0  kcal/ 
mole.  There  has  also  been  the  suggestion  of  a  possible  phase  transition 
in  the  10°-100°K  temperature  range.  The  neutron  diffraction  data  of 
the  present  work  shows  that  NH4CIO4  retains  Its  orthorhombic  (Pnma 
space  group)  structure  from  298°  to  108K.  More  interestingly,  the 
principal  axes  of  rotation  and  the  rotational  amplitudes  of  the  NH£ 

Ions  have  been  extracted  and  are  summarized  in  Table  I. 

Table  I. 

Rotational  Amplitudes  of  NH^  In  NH4CIO4 


Temperature  (8K) 

Princ ipal 
Axis 

Amplitude  (°) 

Direction  Cosines 
a  B  y 

1 

25.7 

.6981  .0000 

.7160 

10 

2 

9.3 

.0000  1.0000 

.0000 

3 

8.6 

-.7160  .0000 

.6981 

1 

30.6 

.6591  .0000 

.7520 

78 

2 

15.8 

.0000  1.0000 

.0000 

3 

15.8 

-.7520  .0000 

.6591 

The  direction  cosines,  a,B,y  are  the  cosines  of  the  angles  between  the 
principal  axes  of  rotation  and  the  x,  y  and  z  crystallographic  axes, 
respectively.  Fig.  1  shows  the  NHj  ion  and  surroundings  at  10°K  as 
viewed  along  principal  axes  1,  the  axis  about  which  the  largest  rota¬ 
tional  amplitude  occurs.  The  78°K  result  is  quite  similar,  with  some¬ 
what  larger  amplitudes.  ’’Tie  ellipsoids  of  thermal  vibration  (Fig.  1) 
indicate  that  H3  and  Hj*  see  identical  potentials  while  Hj  is  more 
weakly  bound  and  H2  is  more  tightly  bound. 

Using  the  model  described  in  Section  II,  QNS  studies  of  NH4CIO3  at 
78°,  105°  and  1308K  have  been  analyzed  in  great  detail  (6).  Below  78°K, 
the  broadening  of  the  QNS  peak  quickly  approaches  that  of  the  instrumental 
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resolution  function.  As  a  first  approximation,  the  QNS  spectrum  at 
each  of  five  Q  values  for  each  temperature  was  least-squares  fit  with 
the  resolution  broadened  theoretical  scattering  cross-section,  adjusting 
Tj  and  T2  (a  T3  -  X3')*  Under  the  assumption  that  the  residence  times, 
T£,  are  related  to  the  rotational  barrier  height  as  x^  ■  x^  exp(+ei/kT), 
activation  energies  ej"  0.2±0.1  kcal/mole  and  C2,3  ■  0.85±0.1  kcal/mole 
are  obtained.  Th*>  bar  denotes  an  average  barrier. 

Raman  scattering  results  (10)  obtained  at  18°K,  49°K  and  82°K 
and  higher  temperatures  for  single  crystal  NH4CIO4  reveal  changes  both 
in  the  NH^  intramolecular  vibrational  modes  as  well  as  the  intermolecular 
(lattice)  modes  generally  consistent  with  the  QNS  and  neutron  diffraction 
results.  Of  particular  interest  is  the  fact  that  three  "soft"  modes 
are  observed  in  this  material. 

Fig.  3  shows  the  spectra  of  the  B^g  and  B3g  symmetry  modes 
which  become  highly  overdamped  as  temperature  is  Increased  from  18°K 
to  82°K.  These  modes  are  Identified  as  NH^  translational  modes  strongly 
coupled  to  rotations.  As  temperature  Increases,  the  rotational  reori¬ 
entations  decrease  the  lifetime  of  the  translational  phonons  leading  to 
the  overdamped  behavior.  Even  more  interesting  is  the  "soft"  mode 
behavior  of  the  Ag  perchlorate  mode  which  occurs  at  45  cm“l  at  298°K. 

The  temperature  dependence  of  this  and  another  Ag  mode  is  shown  in 
Fig.  4.  Based  on  preliminary  neutron  scattering  results  using  the  CNIS 
selection  rules  (2)  the  soft  mode  is  Identified  as  a  perchlorate  trans¬ 
lational  mode.  The  point  to  be  noted  in  this  case  is  that  the  softening 
of  this  mode  begins  to  occur  at  the  temperature  (^100°K)  where  the  NH^ 
reorientation  has  become  very  frequent.  This  mode  continues  to  soften 
as  the  phase  transition  temperature  is  approached.  It  has  been  specu¬ 
lated  that  the  phase  transition  occurs  when  the  perchlorates  begin  to 
"rotate  freely";  the  present  result  is  consistent  with  this  hypothesis. 
Furthermore,  strong  rotation-translation  coupling  is  indicated  for  the 
perchlorate  ions. 

Several  important  conclusions  can  be  drawn  from  our  studies 
concerning  the  molecular  properties  of  NH4CIO4: 

1.  There  is  no  structural  phase  change  from  108K  to  298°K. 

2.  The  observed  dynamics  of  the  four  NH^  hydrogens  are 
consistent  neither  with  all  hindering  potentials  identical  ncr  with  one 
hydrogen  "infinitely"  tightly  bound  as  has  been  assumed  elsewhere. 
Qualitatively,  we  find  e(H^)-c(H3)»E(H3» )-e(H2) ;  quantitatively,  e(H^)- 
0.2+01,  e(H2t3)*0.85±0.1  kcal/mole. 

3.  The  coupling  between  rotational  reorientation  and  trans¬ 
lational  modes  is  manifested  by  the  overdamped  behavior  of  the  NH4  trans¬ 
lational  "soft"  modes;  also  by  the  frequency  shift  of  the  perchlorate 
"soft"  mode.  These  results,  along  with  results  for  azides  discussed  in 
the  next  subsection,  represent  the  first  observation  of  "soft"  modes  in 
energetic  materials. 


PRASK,  CHOI,  IQBAL,  TREVINO 


C.  Metal-azide  Salts 

The  metal-azide  salts  were  the  first  materials  studied  in  a 
systematic  way  within  the  present  program.  (1)  In  terms  of  understanding 
the  Importance  of  anharmonlc  effects  in  energetic  materials,  the  metal- 
azide  salts  provide  numerous  similar  compounds  with  widely  differing 
thermal  and  configurational  stabilities. 

The  main  aspect  of  our  studies  of  metal-azide  salts  that  we 
wish  to  present  here  concerns  the  elucidation  of  anharmonlc  inter¬ 
actions  through  the  study  of  molecular  dynamics  under  the  combined 
effects  of  temperature  and  pressure.  It  should  be  mentioned  that  the 
first  studies  of  this  type  were  performed  only  a  few  years  ago,  and  the 
number  of  compounds  whose  vibrational  properties  have  been  measured 
under  pressure  remains  very  limited.  The  pressure  parameter  is  parti¬ 
cularly  Important  for  the  understanding  of  detonation  phenomena  because 
extremely  high  pressure  precedes  the  propagating  reaction  wave.  Also, 
as  described  in  Section  II,  spectroscopic  studies  performed  as  a  function 
of  temperature  and  pressure  provide  the  only  means  by  which  the  purely 
volume  dependent  anharmonlc  forces  and  the  "self  energy"  can  be  deter¬ 
mined.  Such  studies  have  now  been  performed  for  the  first  time  for  a 
metal-azide  salt,  KN3. 

Fig.  5A  shows  the  pressure  dependence  of  the  Raman-active 
lattice  modes  of  KN3.  (12)  It  is  clear  that  the  rotatory  Eg  and  B^g 
mode  (experimentally  unresolved)  is  much  more  sensitive  to  the  volume 
change  than  is  the  translational  Eg  mode.  Volume  and  single-crystal 
temperature  data  for  the  rotatory  Eg  are  combined  in  Fig.  6.  (12)  In 
this  figure,  VQ  and  o)Q  correspond,  respectively,  to  volume  and  vibra¬ 
tional,  frequency  at  1  atmosphere  pressure  and  298*K;  the  "pressure 
domain"  corresponds  to  increasing  compression  at  constant  temperature, 
the  "temperature  domain"  corresponds  to  Increasing  temperature  at 
constant  pressure.  If  the  self-energy  part  of  the  anharmonlc  inter¬ 
actions  were  small,  the  data  points  would  fall  along  a  single  line  in 
the  temperature  and  pressure  "domains.”  However,  the  data  points  in 
the  temperature  domain  fall  increasingly  far  from  this  line  indicating 
that  the  mode-mode  interactions  involving  the  rotatory  modes  are  quite 
significant.  Aw  is  the  quantitative  measure  of  these  interactions. 

TIN3  is  a  crystal  having  the  KN3  structure  at  298°K.  However, 
with  regard  to  chemical  stability,  TIN3  appears  to  be  transitional  be¬ 
tween  stable  salts  such  as  KN3  and  unstable  salts  such  as  Ag^.  Struc¬ 
turally,  TIN3  is  now  known  to  occur  in  four  different  solid  phases,  all 
reversibly  induced  by  variations  in  temperature  and/or  pressure.  The 
existence  of  the  high-pressure  phase  (III)  was  revealed  through  the 
present  studies.  (13)  It  should  be  mentioned  that  although,  in  general, 
the  correlation  of  physical  and  chemical  instability  has  not  been 
established,  the  results  described  earlier  for  NH4CIO4  suggest  that  at 
least  in  some  cases,  a  strong  correlation  does  exist. 
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In  Fig.  5B  the  pressure  dependence  of  Raman-active  modes  in 
TIN3  up  to  the  II-+TII  phase  transition  (8  kbar)  is  shown.  These  data 
can  be  compared  with  the  results  for  KN3  shown  in  Fig.  5A.  In  the  TIN3 
case,  a  strong  pressure  dependence  for  both  the  rotatory  and  transla¬ 
tional  Eg  symmetry  modes  occurs.  This  is  attributed  to  a  strong 
coupling  between  these  modes  which  reduces  the  R(Eg)  mode  frequency  and 
also  Induces  a  pressure  dependence  in  the  T(Eg)  mode,  compared  to  KN3. 
Since  no  analagous  phase  transition  occurs  in  KN3,  one  must  conclude 
that  the  T(Eg)-R(Eg)  anharmonlc  coupling  in  TIN3  is  a  critical  phase 
stability  parameter. 

With  regard  to  chemical  stability,  an  indication  of  the  differ¬ 
ences  in  interatomic  forces  for  different  azides  is  illustrated  in 
Fig.  7.  The  azide-ion  bending  mode  frequencies  (v2>  are  plotted  vs  the 
heat  of  formation  (AHf)  for  several  metal-azide  salts.  (14)  The  thermal 
stability  of  these  salts  decreases,  approximately,  with  increasing  heat 
of  formation.  The  V2  bending  mode  frequency  clearly  softens  with 
decreasing  thermal  stability  suggesting  a  correlation  between  chemical 
stability  and  mode-mode  Interactions. 

IV.  CONCLUSIONS 

The  present  study  represents  an  important  advance  toward 
achieving  the  goals  of  the  present  program.  In  earlier  work,  structural 
and  molecular  dynamical  properties  were  characterized  only  at  atmos¬ 
pheric  pressure  and  ambient  temperature.  Some  speculation  was  made 
concerning  changes  that  might  be  expected  as  conditions  were  varied  but 
no  concrete  results  were  obtained  for  verification.  In  the  present 
work,  vibrational  properties  of  selected  materials  have  been  measured 
over  a  broad  range  of  temperature  and  pressure.  Rather  significant 
results  are  obtained  including  the  actual  observation,  for  the  first 
time,  of  "soft"  modes  in  energetic  materials.  These  modes  represent  a 
direct  manifestation  of  anharmonlc  mode-mode  coupling,  and  it  is  this 
type  of  coupling  which  provides  the  mechanism  for  intermolecular-  to 
intramolecular-mode  energy  exchange  which  enhances  the  probability  of 
bond  dissociation.  It  is  thus  of  very  significant  interest  for 
energetic  materials. 

Furthermore,  the  present  work  suggests  that  rotational  degrees 
of  freedom  may  play  a  more  important  role  in  the  physical  and  chemical 
stability  of  solids  than  has  been  believed.  This  is  brought  out  in  the 
case  of  NH4NO3  where  the  reorientational  freedom  of  the  ammonium  ions 
is  directly  related  to  the  existence  of  nhase  III.  In  NH4CIO4,  the 
observed  perchlorate  soft  mode  shows  a  softe.iing  correlated  with 
increasing  NH4  rotational  freedom,  and  the  behavior  of  the  mode  appears 
related  to  the  occurrence  of  the  phase  transition.  Such  a  correlation 
is  definitely  established  in  the  case  of  TIN3  where  strong  rotational- 
translational  mode  coupling  is  observed  with  definite  mode  softening  as 
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the  phase  transition  is  approached.  With  regard  to  chemical  instability, 
the  marked  dependence  of  the  rate  of  decomposition  on  the  crystal 
structure  of  NH4CIO4  must  be  related  to  changing  molecular  dynamics. 

Among  the  azides,  the  correlation  of  the  azide-ion  bending  mode  and 
thermal  instability  is  established. 

In  each  of  these  complex  materials  it  appears  likely  that  the 
macroscopic  properties  of  Interest  i.e.  chemical  and  structural 
stability,  are  controlled  primarily  by  a  few,  critical  microscopic 
parameters.  Although  not  yet  explored  systematically,  we  expect  that 
by  properly  choosing  an  additive  or  impurity  atom,  control  of  the 
critical  parameters  could  be  achieved.  This  is  an  area  now  ripe  for 
experimental  investigation.  Furthermore,  recent  results  obtained  by 
computer  molecular  dynamics  studies  (15)  Indicate  that  the  coupling  of 
this  technique  with  state-of-the-art  experimental  work  should,  in  time, 
provide  a  complete  understanding  of  the  properties  of  interest. 

In  summary,  the  present  work  shows  that  the  capability  of 
examining  certain  microscopic  properties  of  energetic  materials  has 
been  advanced  considerably.  Move  importantly,  correlations  between 
molecular  dynamics,  anharmonicity ,  and  millcarily  important  macroscopic 
properties  are  clearly  established.  Finally,  the  comprehensive  approach 
employed  in  these  studies  has  yielded  results  of  significance  not  only 
for  energetic  materials  but  for  all  complex  materials. 
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FIGURE  1:  Nfy*  ion  and  surroundings  in  NH4CIO4  at  10#K  as  viewed 
along  principal  axis  1  (see  text).  A  rotary  displacement  can  occur, 
for  example,  about  the  Hl-N  bond  such  that  H2+H3  site,  H3-+-H3'  site  and 
H3'-*H2  site  without  changing  the  overall  molecular  configuration.  The 
ellipsoids  represent  the  r.m.s.  amplitudes  of  vibration. 


FIGURE  2:  Full  widrh  at 
half  maximum  vs  momentum 
transfer  for  NH4NO3  in 
phases  II  and  IV.  The 
prime  denotes  the  dried 
sample. 
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FIGURE  4:  Temperature 
dependence  of  the  Ag 
symmetry-species  CIO^ 
vibrational  modes.  The 
orthorhombic  to  cubic 
phase  transition  is  at 
513°K.  to,  here,  represents 
the  peak  position. 
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FIGURE  6:  Log  <jj/u)q  versus 
log  v/vD  for  the  rotary  Eg 
mode  In  KN3.  (See  text  *or 
details. ) 
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fjgure  7:  Bending  mode  (Vj) 
O  frequencies  versus  the  heat 

of  formation  for  a  number  of 
metal-azide  salts. 
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